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Finite temperature QCD: progress and outstanding problems 

Carleton DeTar 

Department of Physics, University of Utah 
Salt Lake City, UT 84112, USA 

I review recent progress in numerical simulations of finite temperature quantum chromodynamics and discuss 
the status of some outstanding problems. Included is (1) a discussion of recent results determining the temper- 
ature of the "phase transition" in full QCD, (2) a scaling analysis of the Polyakov loop variable, leading to the 
determination of a constituent quark free energy, (3) studies of critical behavior near the phase transition in two- 
flavor QCD, (4) a discussion of problems and new results in thermodynamic simulations with Wilson fermions, 
(5) recent results in pure gauge theory with a mixed fundamental/adjoint action, and (6) the nonperturbative 
determination of the equation of state with dynamical fermions included. Finally I mention briefly new develop- 
ments in efforts to construct a phenomenology of deconfinement and chiral symmetry restoration, namely (7) the 
dual superconducting model and (8) the instanton model. 



1. INTRODUCTION 

This year has brought new, preliminary re- 
sults from thermodyamic simulations with two 
quark flavors in the staggered fermion scheme at 
Nt = 12, new insights into and questions about 
the critical behavior at the phase transition with 
two quark flavors, a new nonperturbative deter- 
mination of the equation of state for both pure 
Yang-AIills theory and QCD with two sta.ggered 
fermions, and new results for Wilson thermody- 
namics with three quark flavors and at smaller 
lattice spacing, to mention a few highlights. Be- 
cause companion reviews in this volume deal with 
Wilson thermodynamics[l] and the equation of 
state in pure Yang-Mills theory [2], I will treat 
these important topics only briefly. With stag- 
gered fermion thermodynamic data available now 
over a wide range of Ni values it is tempting to 
try to systematize them through scahng relations. 
I offer an attempt toward this end. 
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Figure 1. Polyakov loop susceptibility vs 
from Karsch and Laermann. Solid symbols are 
directly simulated. Open symbols are derived by 
reweighting. 



2. TEMPERATURE OF THE PHASE 
TRANSITION 

Of particular phenomenological interest is the 
temperature of the phase transition or crossover 
from the low to the high temperature regimes in 
QCD. With two flavors of light quarks in the stag- 



gered fermion scheme, there appears to be a dra- 
matic crossover, but so far no evidence for a gen- 
uine phase transition. A phase transition is ex- 
pected at zero quark mass, however[3]. Thus the 
pseudocritical temperature or crossover temper- 
ature is defined as the temperature of the peak 
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Figure 2. Polyakov loop vs temperature. The 
line segments indicate the range of uncertainty in 
locating the maximum slope. 



of a suitable susceptibility, or the temperature of 
maximum change in an observable. Of course, 
the two definitions do not necessarily agree unless 
the chosen susceptibility is rigorously the deriva- 
tive of the observable whose maximum slope is 
sought. Karsch and Laermann use the Polyakov 
loop susceptibility 

XL=N![({RePr)-(RePn (1) 

Combined with a Ferrcnbcrg-Swendsen reweight- 
ing analysis to interpolate between simulation 
points, the method permits a clean determina- 
tion of the pseudocritical couphng, as shown in 
Fig. 1[4]. At larger lattice volumes, however, 
reweighting requires more closely spaced simula- 
tion points and becomes infeasible. The maxi- 
mum slope method is used instead in Fig. 2 to 
locate the phase transition for a wide range of 
lattice sizes[4-7]. Two variables, namely (ReP) 
and (ijjijj) were examined to determine the max- 
imum slope. This figure shows results only for 
the lowest available quark mass in each data set. 
Included are new results at A'"^ = 8 and pre- 
liminary results at Nt = 12 [7], the latter corre- 
sponding to a lattice spacing of 0.1 fm. The tem- 



perature scale, based on the rho meson mass, is 
discussed in greater detail below. The tempera- 
tures thus determined are summarized in Fig. 3. 
There is an encouraging consistency in the re- 
sults, which place the crossover temperature at 
approximately 140-160 MeV, based on the p me- 
son mass. It should be emphasized, however, that 
over the range of lattice parameters of this com- 
pilation, the nucleon to rho mass ratio is approx- 
imately 1.5, 20% above the physical value. Us- 
ing the nucleon mass to set the temperature scale 
would therefore result in a 20% reduction in the 
crossover temperature. Physical values of this ra- 
tio have been reported in quenched calculations 
for both Wilson and staggered fermions[8], offer- 
ing hope for future simulations with dynamical 
fermions. Furthermore, the mass ratio 111^^ /nip 
lies in the range 1/3 to 1/2, twice to three times 
the experimental value, because the quark mass 
is too high. Indeed, an extrapolation to zero 
quark mass was attempted in constructing Fig. 3. 
However, present indications suggest that lower- 
ing the quark mass does not substantially change 
the crossover temperature. 

The temperature scale in Fig. 2 was con- 
structed following Blum et al [6,9] and is based 
on the lattice rho meson mass. That is, from 
several zero temperature simulations[12] an in- 
terpolation fitting formula is constructed, giving 
amp{6/g^ ,'niqa). For present purposes, we use 
the combination of cubic splines given by 

amp = So{6/g'^) + Si{6/g'^)amg 

+S2{6/g')(am^,r (2) 
m^/nip = (am<,)^/25o(6/g^) 

+Si{6/g'^)amg (3) 
+S2i6/g'){amgf^' 

with knots at G/g"^ = 5.3, 5.5, and 5.7. Fig- 
ure 4 shows the parameter range of thermo- 
dynamic simulations and companion zero tem- 
perature spectrum simulations used in this re- 
view, indicated by the plot symbol "s". As is 
usual, the thermodynamic simulations outpace 
the spectrum simulations. Thus for 6/g^ > 5.7 
encountered in the Nt = 12 simulations, val- 
ues higher than any available zero temperature 
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Figure 3. Temperature of the crossover in units 
of the rho meson mass with two Hght quarks vs 
the squared ratio of the pion to rho mass. For 
points from the staggered simulations, the mass 
of the second (non-Goldstone) pion is used. The 
curved line segment shows the error bar for the 
new Nt = 8 staggered point. The vertical dashed 
line indicates the physical mass ratio. 



spectrum study, an extrapolation becomes nec- 
essary. For this purpose I used an approximate 
tadpole-improved asymptotic scaling formula[10], 
based on high temperature plaquette measure- 
ments. A temperature in MeV is then inferred 
from T = 770/[amp{6/g'^,mga)Nt]. This scale 
does not allow for variations in the rho mass with 
quark mass. It also inherits all the systematic er- 
rors of the zero temperature spectral calculations. 
Nonetheless, the rough figure 140-160 MeV for 
the crossover is already extremely useful for mod- 
els of quark plasma formation. 

3. CONSTITUENT QUARK FREE EN- 
ERGY 

The wide range of Nt values available this 
year makes possible an amusing analysis of 
the Polyakov loop variable, which measures the 
change in the free energy of the thermal ensem- 
ble due to the introduction of a point spinless test 
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Figure 4. Thermodynamic simulation points 
(conventional plot symbols) and spectrum simu- 
lation points (letter s) for staggered fermion data 
used in this review. 



quark. This free energy difference 
/(T,m,) = -Tlog(ReP/3), 



(4) 



a function of the temperature T and light 
quark mass m,, includes the lattice- regulated 
ultraviolet-divergent self-energy of the point 
source, proportional to the inverse lattice spacing 
1/a, and the free energy of the screening cloud 
of light antiquarks and quarks, which we might 
call with some risk of confusion the "constituent 
quark free energy" . Computing the self energy to 
leading order in perturbation theory, we have 



/(T, m,) = 2wCFav7/a + /cq(T, m,). 



(5) 



where Cp = 4/3 is the color Casimir factor for 
the triplet representation, ay is the color fine 
structure constant for appropriate to heavy quark 
bound states at the same lattice scale, and 
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A^f ^ 6 - 2 cos(27rfc^/Ar,) + {Maf 



(6) 



is the dimensionless static lattice propagator for a 
Debye-screened electrostatic gluon field evaluated 
at zero separation. Why screening? Although the 
ultraviolet divergent contribution 1/a is imiquely 
determined in the limit a — » 0, an arbitrary in- 
frared cutoff, here embodied in the Debye mass 
M, determines where the contribution from the 
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point quark ends and the contribution from the 
screening cloud begins. Thus there is no unique 
definition of "constituent quark free energy" . In- 
stead, within the framework of any consistent def- 
inition, Eq (5) permits a separation of two contri- 
butions, one varying in a known way with the lat- 
tice scale 1/a and the other, unknown, but scale- 
invariant. Herein Hes its predictive power. 

For present purposes, using Nt = i/aT, I chose 
a simpler approximate form 



h^{T,m,) = -T(log(ReP/3) + cNt) 



(7) 



and adjusted the dimcnsionlcss constant c by eye 
to achieve the rough scaling agreement shown in 
Fig. 5. For each data set only values for the light- 
est available quark mass are used. Although the 
quark mass values niq/T are not the same from 
one Nt to the next in this figure, they are small 
{'niq/T < 0.1) and would be expected to con- 
tribute little (of the order 10 MeV) to the free 
energy. Thus one would expect only a small in- 
consistency from this variation. The best value 
for c appears to be about 0.4, compared with 
values ranging from 0.26 to 0.35 expected from 
the lowest order perturbative self energy with a 
screening mass M = 3.2T. We note that at the 
crossover, the free energy drops by about the 300 
MeV expected in a constituent quark model with 
deconfinement at high temperature. 

4. CRITICAL BEHAVIOR 

The phase structure of QCD also has impor- 
tant phenomenological consequences. Sketched 
in Fig. 6 is a proposed phase diagram with two 
flavors of light quarks and one strange quark 
(2-1-1 flavors), showing in what mass range a 
thermal phase transition is expected, and whether 
it is flrst or second order. Such a phase struc- 
ture is motivated by an analysis of the cor- 
responding sigma models in mean fleld theory, 
augmented by an analysis of quantum fluctua- 
tions[3,ll]. Whether QCD conforms to this ex- 
pectation remains to be established. Simulations 
with 2 + 1 flavors in the staggered and Wilson 
fermion schemes both support the existence of 
the first order region, but do not agree on its ex- 
tent[l,14,15]. Staggered fermion simulations of 
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Figure 5. Constituent quark free energy as de- 
fined in text. 



the Columbia group found evidence for a first or- 
der signal for (m„,d,ms) « (15,15) MeV, but 
no signal for {mu,d,'ms) « (15,30) MeV. A re- 
cent Wflson fermion simulation by the Tsukuba 
group found a first order signal with quark masses 
as large as {mu.d.fns) « (150,150) MeV and 
(m„,d,m,) « (0, 400) MeV[l,15]. Is the discrep- 
ancy a strong coupling artifact? At strong cou- 
pling the fermion doublers in the Wilson scheme 
tend to increase the effective number of relatively 
light flavors. With more flavors the chiral phase 
transition is stronger. On the other hand in the 
staggered fermion scheme, the breaking of fla- 
vor symmetry in strong coupling may reduce the 
efl^"ective flavor number, thereby weakening the 
phase transition. 

A recent analysis of the three-dimensional 
Gross-Neveu model by Kocic and Kogut ques- 
tions the conventional wisdom that places QCD 
with its composite scalar mesons in the same uni- 
versality class as sigma models with their ele- 
mentary scalar mesons [16]. One would expect 
that the 3D Gross-Neveu model would exhibit 2D 
Ising universality. Instead, in a detailed simula- 
tion Kocic and Kogut found mean-fleld scaling. 
Now for QCD and the quark plasma there may be 
little phenomenological importance to such sub- 
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Figure 6. QCD phase diagram for 2 + 1 flavors 
as a function of a degenerate up and down quark 

mass m^u.d = n^u = fnd and a strange quark mass 
Wis. The heavy line is the critical phase boundary. 
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Figure 7. Crossover coupling vs quark mass for 
two flavors of staggered quarks from Karsch and 
Laermann. Solid Hne is a fit with 0(4) critical 
exponents. Dashed line is a free fit. Diamonds 
are new results. 



tleties, but the question is of broad significance 
for understanding critical behavior in field theo- 
ries. Furthermore, to take the continuum limit of 
lattice simulations, we must learn to do extrap- 
olations to small quark mass and large volume. 
In the vicinity of the phase transition, a correct 
extrapolation requires the correct critical expo- 
nents. Thus an important test of this proposed 
phase structure is to determine whether QCD has 
the critical behavior expected from universality 
along the critical phase boundary. 

For the remainder of this section I will focus 
on the two flavor theory in the staggered fermion 
scheme, corresponding to the upper portion of the 
TRs axis in Fig. 6. In this case 0(4) universality is 
expected in the continuum limit. To be more pre- 
cise, in the staggered fermion scheme, one expects 
0(2) critical behavior on coarse lattices where the 
flavor symmetry breaking of the staggered scheme 
is significant, and 0(4) as the continuum limit is 
approached. Recent work by Karsch and Karsch 
and Laermann attempts a determination of some 
of the critical exponents of QCD[4,17]. They ex- 
ploit the standard sigma model analogy between 
QCD and a magnetic system. In this analogy the 



quark mass plays the role of an external magnetic 
field and (^^) plays the role of the magnetiza- 
tion. 

The standard analysis of critical behavior be- 
gins with the assumed scaling of the critical con- 
tribution to the free energy in a magnetic sys- 
tem[18]. This contribution is singular and dom- 
inant at small field close to the zero field criti- 
cal temperature Tc(0). The scaling property, ex- 
pressed in terms of the scaled temperature t = 
[T - T,(0)]/T,(0) and magnetic field h, is 

/cnt(i,/l) =fo-Vcnt(fo^'i,fo^'-/l). (8) 

From this scahng behavior one derives a scaling 
relation for the critical contribution to the mag- 
netization s = —dfciit/dh: 

s{t,h) = h'l'y{x) (9) 

where x = th^l^^ and y{x) is a scaling function 
(often called an equation of state in the statisti- 
cal mechanics literature). Here d and 6 are stan- 
dard critical exponents derived from yt and yu- 
An important consequence of this result is that a 
crossover peak in the susceptibility Xh = ds/dh 
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moves along a curve of constant x = Xpc as h and 
t are varied. Thus if critical scaling holds, once 
the pseudocritical temperature is known at one 
h, it can be predicted at all h. 

In QCD the quark mass plays the role of 
the magnetic field and {il>ip), the magnetiza- 
tion. Specifically, Karsch suggests using h = 
rUg/T = arUgNt and i = 6/g^-6/g'^{0,Nt), where 
gc{0,Nt) is the critical gauge coupling at zero 
quark mass for a particular A''([17]. Scaling then 
implies that 



6/gL{mga) = 6/gl{0) + {mg/T) 



(10) 



Using this expression Karsch presented an anal- 
ysis of the crossover for Nt = 4, 6, and 8 for two 
quark flavors[17]. He found encouraging agree- 
ment. The addition of new data this year at 
Nt = A permits a more refined analysis, shown 
in Figure 7[4]. Their best fit critical exponent 
l/PS is 0.77 ± 0.14, in slight disagreement with 
the 0(4) value 0.55(2), but consistent with the 
0(2) value 0.60(1) and the mean field value 0.67. 

Karsch and Laermann also introduced a new 
cumulant 



A 



d\n{ipip{6/ g'^ ,mq)) 
dlnvia 



xy'{x) 



(11) 



S 36y{x) 

that evaluates the critical exponent A = 1/^ at 
X = t = 0. They obtain 0.21 < 1/^ < 0.26 
consistent with the 0(4) value 0.208(2) and 0(2) 
value 0.2080(3), and somewhat inconsistent with 
the mean field value 0.33. Thus the Kocic-Kogut 
scenario cannot be decisively excluded. 

To carry these results further, I test the scal- 
ing relation (9) in simulations with two flavors 
of staggered fcrmions over the wide range of cur- 
rently available Nt. For this purpose I use shghtly 
diff'erent variables to permit comparison among 
diff"erent Ni and to avoid quantities with anoma- 
lous dimensions, namely, 

h = rnl{mg,T = 0)/ml{m,T = 0) (12) 
t = [T-T,(0)]/T,(0) (13) 
s = h-'^m{i>i>{mg,T))/T^ (14) 

The scaling relation (9) then gives a universal 
function 
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Figure 8. Scaled (V'V') '"^ scaled temperature with 
0(4) critical exponents. 



with X = th^l^^ . The extra factor h^^ is needed 
to compensate the quark mass factor nig. Let us 
apply this analysis to data for (i/'^) from several 
groups[5]. Setting the critical exponents 6 and 
(3 to their 0(4) values and adjusting the sole re- 
maining free parameter Tc(0) in 10 MeV incre- 
ments to get the best agreement I get the re- 
sult shown in Fig. 8. With the exception of the 
A''^ = 12 data, the scaling agreement is rather 
good. At this level it is not possible to distin- 
guish 0(4) from 0(2) and mean fleld critical be- 
havior. Setting Tc(0) to 140 MeV or 170 MeV 
worsens the agreement noticeably, but 160 MeV 
gives comparable consistency. Obviously a host of 
systematic errors, including finite volume efli"ects 
and deviations from continuum scaling enter the 
analysis, so refinements are certainly needed be- 
fore the method can serve as a definitive test of 
critical behavior. The most glaring inconsistency 
in this figure comes from the preliminary Nt = 12 
data. Increasing TdO) for only this data set to 
160 MeV and plotting it with other data com- 
puted for Tc(0) = 150 MeV brings the Nt = 12 
data at the lower quark mass into good agree- 
ment. Thus the discrepancy could be caused ei- 
ther by a gradual upward shift in the crossover 
temperature as the lattice spacing is decreased. 
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or by an erroneous extrapolation of the hadron 
spectrum above = 5.7, or by a failure of the 
scahng hypothesis over this parameter range. 

5. THERMODYNAMICS WITH WIL- 
SON FERMIONS 

In order to claim we can infer features of contin- 
uum QCD from numerical simulations it is essen- 
tial that we demonstrate that our answers are in- 
dependent of the fermion scheme. Unfortunately, 
thermodynamic simulations with Wilson fermions 
are not sufficiently developed at present to make a 
confident comparison with the staggered scheme. 
The fundamental difficulty is that we are deal- 
ing with a phase transition associated with the 
restoration of a spontaneously broken chiral sym- 
metry, but the Wilson scheme builds in an explicit 
breaking of this symmetry, which goes away only 
in the continuum limit. 

It is popular to define a chiral lino k = Kc(6/g^) 
where a number of characteristics of chiral sym- 
metry (vanishing zero temperature pion mass, 
vanishing current quark mass) hold, at least with 
some consistency. Thus we are interested in 
studying the crossover or phase transition in a re- 
gion close to the chiral line. The Tsukuba group 
found that to do so requires either working at 
very small values of Q/g"^, usually considered to 
be in the strong coupling regime, or working at 
much higher Nt than is feasible with present re- 
sources[19]. The former choice risks encountering 
lattice artifacts, and the latter is expensive. In- 
deed, evidence for lattice artifacts in two flavor 
simulations at A^t = 6 was recently reported by 
the MILC collaboration, which found a first or- 
der, possibly bulk, phase transition in simulations 
at K = 0.17, 0.18, and 0.19 in close proximity to 
the thermal crossover [20]. The Tsukuba group 
has continued its pioneering work at small 
For an authoritative review, see Iwasaki's talk[l]. 

At this conference results of the first simula- 
tion at A^t = 8 in the Wilson scheme were re- 
ported by the MILC collaboration [21]. The sim- 
ulation was done at = 5.3 over a range of 
K up to Kc « 0.168, allowing a comparison with 
previous results at A^t = 4 and 6. The thermal 
crossover, now shifted to Kt ~ 0.167, shows no ev- 
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Figure 9. Chiral condensate vs k for two fiavors 
of Wilson fermions at Q/g^ = 5.30. 



idence for the metastability seen at slightly larger 
K at Nt = 6. As illustrated for {ijjil') in Fig. 9, as 
Nt is increased, bulk quantities appear to follow 
an envelope established by the zero temperature 
theory, breaking away at the crossover. More- 
over {'ip'tp) appears to be decreasing immediately 
prior to the crossover, suggesting progress toward 
a low quark mass. For details, please see Tous- 
saint[21]. Thus, fortunately, the disease seen at 
slightly larger k in the iV^ = 6 simulations appears 
not to be spreading to lower k as Nt increases. 

It is clear that with present methods, Wilson 
fermion thermodynamics are far more costly than 
staggered fermion thermodynamics, requiring far 
larger lattices to suppress strong-couphng arti- 
facts. Further progress with the Wilson scheme is 
likely to require working with an improved action. 
Indeed the Tsukuba group has adopted one such 
improvement and they report some preliminary 
results at this conference[l,22]. 



6. PURE YANG MILLS WITH MIXED 
ACTION 

This year Gavai, Grady, and Mathur reported 
results for simulations of SU(2) Yang-Mills theory 
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with a mixed fundamental/adjoint action 

S = 0fJ2[^-j^ReTifUp] 

p 

+0aJ2il-^TrfUlTvfUp]. (16) 
p 

where = 2 for SU(2), /?/ and /3a are the fun- 
damental and adjoint gauge couplings, the gauge 
link matrices U are in the fundamental represen- 
tation, and Tr/ is the trace in the fundamental 
representation [23]. They have reopened the ques- 
tion of the interplay between the well-known bulk 
transition[25] and the thermal phase transition. 

Continuum universality requires that as Nt is 
increased, the entire thermal phase boundary sep- 
arating the confined and deconfined phases must 
shift toward weak coupling. If instead one end 
remains anchored to a bulk phase boundary, the 
location of which by definition is asymptotically 
constant as A^t — » oo, it would be possible to 
approach the continuum limit at zero tempera- 
ture along some directions in the Pa — /3/ plane, 
all the while remaining in the deconfined phase, 
and to approach the continuum limit at zero tem- 
perature along other directions in the confined 
phase. Gavai, Grady and Mathur found that at 
Nt = 4: the second order thermal phase boundary 
in SU(2) connects smoothly to what had been 
hitherto called the first order bulk phase bound- 
ary. If universality were to survive, the Nt = 6 
phase boundary would have to show a shift to- 
ward weak coupling. Indeed, in a more recent re- 
markable work Mathur and Gavai found that at 
Nt = 6, the strangely coalesced phase boundaries 
appear to have shifted together toward weak cou- 
pling slightly [26], suggesting, perhaps, that the 
first order segment had been mislabeled as a bulk 
phase transition. 

At this conference Heller presented results of 
recent simulations in SU(3) Yang-Mills theory, for 
which an apparently different and perhaps less 
surprising picture is emerging[27]. Results are 
summarized in Fig. 10. For Nt = 4 the thermal 
and bulk phase transitions are found to coalesce 
in SU(3) as in SU(2). But for TV* = 6 and 8 
the thermal phase boundary clearly breaks away 
from the terminal segment of the bulk line, joining 
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Figure 10. Phase diagram for pure Yang-Mills 
theory with a mixed fundamental-adjoint action. 
The bulk phase boundary, applicable to A^j = 4, 
6, and 8, is plotted as a darker line. The burst 
marks the bulk endpoint. The Nt = 8 points were 
found from an infinite volume extrapolation. 



it above the bulk endpoint. As Nt is increased, 
the junction moves to higher ffa, so that the en- 
tire thermal line shifts toward weak coupling in a 
manner consistent with universality. It would be 
surprising if the SU(2) result were to remain so 
dramatically different at still higher Nt. 

The signature for a separation between the bulk 
and thermal phase transitions in mixed action 
Yang-Mills theory is similar to what was found 
with Wilson fermions at large k in the A^t = 6 sim- 
ulations, namely, with increasing df, a first order 
phase transition is first encountered, signaled by 
a discontinuity in the plaquette, but without an 
appropriate jump in the Polyakov loop variable. 
At slightly higher /3f the Polyakov loop variable 
then rises, signahng a thermal phase transition or 
crossover. Now the Wilson fcrmion determinant 
surely induces an adjoint plaquette, which might 
become strong at large k. Thus there may be 
more than coincidence relating the two plienoni- 
ena. However, at this conference Runiniukainen 
reported simulations measuring the strength of 
induced fundamental and adjoint couplings using 
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a microcanonical demon method. He found a dis- 
appointingly small induced adjoint action, raising 

doubts about such an explanation for the pecu- 
liar behavior of Wilson thermodynamics at large 

k[28]. 

7. EQUATION OF STATE WITH NON- 
PERTURBATIVE METHODS 

The equation of state for hadronic matter, 
giving the energy density e(T, m,) and pressure 
p{T,mq) as a function of temperature and quark 
mass is another quantity of phenomenological im- 
portance. The earhest determinations of these 
quantities used the basic thermodynamic identi- 
ties 

= am ^''^ 

dF 

av 

In a lattice simulation each such derivative of the 
free energy involves a separate variation of the 
spatial and temporal lattice constants as and at, 
requiring a concomitant renormalization of the 
gauge coupling. Some years ago Karsch deter- 
mined the perturbative asymptotic variation of 
the gauge coupling with respect to the anisotropy 
parameter ^ = ai/as[29]. Unfortunately, present 
simulations are not in the perturbative asymp- 
totic region. Although in principal the Karsch co- 
efficients could be determined nonpcrturbativcly 
from lattice simulations, this has not yet been 
done successfully [9, 30]. 

Fortunately there is a different nonperturbative 
route to e and p. The so-called interaction mea- 
sure 



PIT 



(18) 



I = e-2,p 



(19) 



is more easily determined, since it involves an 
isotropic variation of the lattice constant, requir- 
ing only the usual nonperturbative renormaliza- 
tion of the lattice coupling. The pressure, on the 
other hand, can be determined separately from 
the free energy / = —p by integrating either of 
two relations[31] 



1 



dF 



1 dF 
V dm„ 



(21) 



A vacuum subtraction is performed to give the 
pressure relative to the pressure of the nonper- 
turbative vacuum: 



P_ 

rp4 



P_ 
rp4 



27V/ / ' rf(6/g'^)[(□(6/g'^am,)) 

./cold '■ 

- (□(6/5'2,am,))J (22) 

./cold ^ 



{ijjijj{Qlg'^,am'g)) 



(23) 



21/ (9(6/g2) 



(20) 



Of course the latter equation may be used only 
when dynamical quarks arc present. Together 
with a determination of the interaction measure, 
this result can then be used to determine the en- 
ergy density. 

This year new nonperturbative determinations 
were reported for the energy density and pressure 
in SU(2) and SU(3) Yang-Mills theory[2,32] and 
in SU(3) with two flavors of dynamical staggered 
fermions[6,33]. In the latter case, although there 
is no evidence for a bona fide phase transition at 
nonzero quark mass, nonetheless there is a steep 
rise in the energy density at the temperature asso- 
ciated with the largest slope in the Polyakov loop 
and (V'V')i seen in Fig. 11. In a coohng quark 
plasma such a strong crossover could cause a mo- 
mentary slowing in the expansion of the plasma 
as the quarks and gluons reorganized themselves 
into more compact hadrons. 

8. MAGNETIC MONOPOLES AND DE- 
CONFINEMENT 

A study of the thermal behavior of QCD may 
provide insight into the mechanism of confine- 
ment. An old 't Hooft — Mandelstam model char- 
acterizes the confining QCD vacuum as a dual su- 
perconductor, with an electric Meissner effect and 
a condensate of color magnetic monopoles[34]. 
Some years ago Schierholz et al and Kronfeld 
et al explored the association between confine- 
ment and the presence of monopole currents in 
Yang-Mills theories[35]. Interest has revived re- 
cently [36, 37]. 
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Figure 11. Energy density (upper curve) and 
three times the pressure (lower curve) vs tem- 
perature (scale based on the rho meson mass) 
for two light quark flavors {niq/T = 0.1) in the 
staggered fermion scheme from Blum, Gottlieb, 
Karkkainen, and Toussaint. The bursts give the 
pressure extrapolated to zero quark mass 



To identify monopole currents in a nonabelian 
gauge theory it is necessary to carry out a U(l) 
projection of the gauge links. This is done by 
first fixing a suitable gauge. Popular choices in- 
clude maximal Abelian gauge and a variety of 
"unitary" gauges, one of which involves diago- 
nalizing the product of gauge links forming the 
Polyakov loop. A "monopole" contribution is 
then extracted from the resulting U(l) gauge field 
following the procedure of DeGrand and Tous- 
saint [38] . 

In the new approach the string tension and 
other confinement features are computed using 
only the monopole contribution to the U(l) field. 
Good agreement is found with the full string ten- 
sion computed in the conventional way. In the 
past year the Kanazawa group has also calculated 
the projected U(l)-monopole Polyakov-loop ex- 
pectation value in both SU(2) and SU(3) Yang- 
Mills theory [39]. The behavior of the monopole- 
projected Polyakov loop variable imitates the be- 
havior of the conventional Polyakov loop variable. 
The similarity is seen in a variety of U(l) projec- 
tion schemes. 

Although the results are promising, further 
work is needed (1) to demonstrate that the suit- 



ably defined abelian monopole currents survive 
the continuum nmit[40] and (2) to find a suitable 
order parameter for monopole condensation [41]. 

9. INSTANTON MOLECULES AND 
THE CHIRAL PHASE TRANSITION 

The Stony Brook instanton program classifies 
the dominant gauge field configurations in the 
thermal ensemble in terms of their instanton con- 
tent and seeks to understand the principal fea- 
tures of QCD in terms of them. In recent work 
Ilgenfritz and Shuryak and Schafer, Shuryak, and 
Verbaarschot argue that in full QCD, as the chiral 
phase transition is approached, the light fermion 
determinant induces an attractive interaction be- 
tween instantons and anti-instantons[42]. The re- 
sulting molecules are predicted to predominate 
over soHtary instantons and anti-instantons. It 
is argued that the strong correlation drives the 
rliiral phase transition. Schafer, Shuryak, and 
Verbaarschot compute hadronic screening masses 
in the model and find a spectrum in qualitative 
agreement with results from lattice simulations. 
It would be interesting to test their proposals fur- 
ther in lattice simulations. 
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